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ABSTRACT 


A ten-level primitive equation numerical model is used to study the influences of 
the bottom topography in the California Current System. Five different numerical 
experiments were integrated for 40 days after being initialized with a baroclinic along- 
shore coastal jet representative of the observed coastal jet. By comparison with a flat 
bottom case, in Experiments | and 3, in which the topographic slopes face westward or 
northward, the topography appears to have stabilizing influences on the mean flow. In 
Experiments 2 and 4, in which the topographic slopes face southward, the topographic 
§-effect appears to be strong enough so that these two cases have already reached a 
quasi-steady state by 40 days. This is because in Experiments 2 and 4 the long, non- 
dispersive topographic Rossby waves are effective at transporting the eddy energy away 
from the source region near the coast. The resulting new mean flow appears to be 
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I. INTRODUCTION 


The eastern boundary current regime off the coast of California is characterized 
by the development and evolution of eddy activity and by a continual interaction 
between this eddy activity, the “mean” California Current System (CCS), and important 
physical processes in the mixed layer. Baroclinic instability of the seasonal mean flow 
has been identified as an important generating mechanism in the Northeast Pacific 
Ocean, and off Vancouver Island (Wright, 1980; Emery and Mysak, 1980; Thomson, 
1984). It has also been shown to be important off the coast of Northern California and 
Orgeon (Ikeda and Emery, 1984), and it may also be an important mechanism in the 
central California coast region as well. The pronounced mesoscale eddy activity in the 
CCS has been identified and best seen in satellite IR images, especially those from the 
AVHRR on board NOAA-6. An eddy seen in such satellite data during the CODE 
experiment off Vancouver Island (Thomson, 1984) was shown to develop due to 
baroclinic instability and it was observed (in CTD and XBT data) to extend form the 
surface to a depth of 1 km. In the CCS however, the vertical structure of the eddies 1s 
not Well known. Furthermore, the kinematics, the energetics and the vertical structure 
of an eddy may change with time as the eddy evolves and migrates off-shore and 
downstream into the “mean” CCS. 

In most studies of flow over topography, it has been found that the slope effect 
(topographic-B) acts to stabilize the flow. However there are two recent studies of flow 
over topography which indicate that enhanced instability can result. The first study is 
the analytic work of Charney and Flierl (1981), who demonstrated that the stability of 
a parallel mean flow in the presence of topography can be quite different from that of 
the identical mean flow in a flat bottomed ocean. They examined the specific problem 
of destabilization of the mean flow by form drag instability of a baroclinic flow which 
is neutrally stable in the absence of topography. The instability occurs for topographic 
scales of the order of 70-100 km, with growth rates proportional to the topographic 
height. This means that the available potential energy of the mean flow can be trapped 
by the orographic instability even in situations where normal baroclinic instability is 
unable to extract energy from the mean flow. The second study Is laboratory modeling 


by Narimousa and Maxworthy (1985). In their experiments, coastal upwelling is 


produced by an along-shore stress in a rotating cylindrical tank in the presence of a 
cone-shape ridge in the bottom topography. As the upwelling process continues, the 
interface (front) between the two layers of the two-layer fluid used in the tank surfaces 
and moves away from the wall (coast). When the front has reached a certain distance 
from the wall, baroclinic waves develop on the front, and large standing waves develop 
downstream of the topographic ridge. Under some circumstances, cyclones pinch off, 
both from the upstream front at the topographic ridge and from the crests of the large 
standing waves downstream of the ridge. These laboratory results show many features 
reminiscent of eddy flows recently observed in the CCS. The results of these 
investigations suggest that a study of the baroclinic and barotropic instability processes 
over a pronounced escarpment normal to the coast is needed to fully understand the 
oceanic eddy field off the coast of central California. 

This study is designed to identity the role of topography in modifying the 
stability of the CCS. An objective is to determine whether the eddy field in the CCS is 
the result of an instability of the mean flow or not, and whether the topography effects 
the instability. This is done by performing initial value numerical experiments with an 
existing multi-level primitive equation model adapted to the CCS. The model domain 
covers 6° longitude by 6° latitude and has approximately 8 km X 10 km horizontal 
resolution with 10 levels in vertical. The space staggered B-scheme in horizontal and a 
6 coordinate system in vertical were used. Five different numerical experiments with 
four different bottom topographies and one flat bottom case were performed. We note 
that the results of the experiments were different and in the results have different 


interpretations. 


Il. MODEL AND EXPERIMENTS DESIGN 


A. MODEL DESCRIPTION 

Te purpose of this study is to investigate the influences of the topography in the 
CCS. The present numerical experiment uses a ten level primitive equation model 
(Haney, 1985). The model is in sigma coordinates (non-dimentional depth) and has 
open boundaries (Camerlengo and O’Brien, 1980) on all but the eastern coastal 
boundary. 

The present domain of the model ocean 1s the rectangular region extending from 
approximately 124° W to 130° W and from 36.5°¢ N to 42.5° N (Fig 2.1). The present 
model has 65 X 65 = 4225 grid points uniformly spaced in the horizontal with a 
resulting grid spacing of 8 km in the east-west direction and 10 km in the north-south 
direction. Figs. 2.2-2.6 show 3-dimensional perspective plots of the actual bottom 
topography of the region and those used in the model ocean. Four numerical 
experiments, each with a different bottom topography (Figs. 2.3-2.6), were run in order 
to study the role an idealized escarpment (Mendoceno Escarpment) on the CCS. 

This study only focuses on the influences of variable bottom topography, because 
a variable bottom topography may be a more important factor than an irrgular 
coastline (Peffley and O’Brien, 1976; Preller and O’Brien, 1980). This conclusion is also 
supported by the recent and highly relevant laboratary modeling work by Narimousa 
and Maxworthy (1985). In the four experiments there was no wind or thermal forcing, 
but the initial temperature field was specified so as to produce a geostophically 
balanced, meridional (long-shore) shear flow having a gaussian-shaped jet. The model 
is based on the primitive equations which are integrated in time starting from the 


specified initial condition. 
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Figure 2.1 Model region. 
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B. EXPERIMENTS DESIGN 
The numerical model was used in this study to carry out four different 
topographic experiments: 
1. off-shore slope only 
2. escarpment with shallowest depths to the north 
3. escarpment with deepest depths to the north 
4. ridge 

Of these configurations, Fig. 2.3 corresponds to the continental rise off 
California, Fig. 2.4 corresponds to an idealized Mendoceno Escarpment, while Fig. 2.5 
and Fig. 2.6 are for comparison with Fig. 2.3 and Fig. 2.4. The laboratory experiments 
of Narimousa and Maxworthy (1985) used a ridge (Fig. 2.6), but an escarpment (Fig. 
2.4) is clearly more realistic of the actual topography in the CCS. 

The initial temperature field approximately followed the Orlanski and Cox (1973) 
paper, which was a study of the Gulf Stream, but modified to represent the CCS. In 
order to define a baroclinic jet corresponding to the CCS, and to fit the OPTOMA 
(Whittman, et al, 1985) profile of buoyancy frequency N? (largest N? at about depth 
100m) and temperature in the upper ocean, we used four analytic expressions for 
temperature. In addition, we disturbed the equilibrium state with a small random 
pertubation in the temperature field. The random number (R) was between -0.05 and 
+0.05 °C, and it was made to decrease in magnitude with increasing depth. We need 
larger surface temperature gradients because there is no wind forcing in this model. The 
realistic surface temperature is 10°C in-shore and 18°C off-shore in summer, so we 
have to adjust T, and T., to get both a realistic surface temperature and geostrophic 
transport as in Freitag and Halpern (1981). In addition, we use T, and T, to fit the 
N? profile in the OPTOMA data. So the analytic expression for the initial temperature 


is defined as follows: 


elie ot Te ee ~ eXp(-zzz/ 450, ) (eqn 2-1) 
Ty =a ele = 2.0C 
ess ex Ljexp(-ex) X DT 


— 
I 


exp (-Zzz/ 150.) 


= 
ll 


AMPI X (-zzz/L + 1 + 1/2n(sin(2m(1/2-zzz/L)))), for zzz less than L 
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Figure 2.2 Actual bottom topography. 
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Figure 2.3  Off-shore slope only. 
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Figure 2.4 Escarpment with shallowest depths to the north. 
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Figure 2.5 Escarpment with deepest depths to the north. 
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eee VIE 2 exp (-222/450.) 

In these expressions x is distance off-shore, zzz is depth (in meters, positive 
downward) and 1/€ is the off-shore length scale for the baroclinic jet (40 km). The 
vertical scale for the jet is 150 m and the vertical scale for the main thermocline (and 
the random disturbance) is 450m. T, is defined so as to produce a realistic alongshore 
geostrophic baroclinic jet,V(x,z). This is done using the hydrostatic and geostrophic 
(thermal wind) relation and the assumption that the vertical average current is zero. 


The constants used in this model are in Table 1. 


TABLE 1 
VALUES OF CONSTANT USED IN THE MODEL 


DTAU = 800. sec integration time step 


G = 981.0 cm/sec? acceleration of gravity 


C,, = 0.958 cal/(gm °K) specific heat of sea water 


We = 273.2 °K reference temperature 

@= 2.0 X 1.0E-4 (1/°k) thermal expansion coefficient 
F = 1.0E-4 (1/sec) coriolis parameter 

Q = 2n day earth rotation rate 

p = 1.0276 gm/cm? 


Doo = 4500m deepest ocean depth 


density 


L = 200m wavelength for the seasonal thermocline 
DT = 11°C horizontal temperature variation 
Ap ly 3. °C amplitude of seasonal thermocline 


Avie, = 12°C amplitude of main thermocline 





Since the slopes of the bottom topographies used in the experiments (Figs. 
2.3-2.6) are not excessive, the resulting model flows can be interpreted (at least 
qualitatively) using quasi-geostrophic theory. From the quasi-geostophic potential 
vorticity equation we can define two parameters y and 6 which describe the potential 


for mean flow instibilities (Killworth, 1980). The parameter y is the ratio of the 
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horizontal length scale of the jet to the Rossby radius of deformation, and 6 is the 
ratio of layer depths (or the fraction of the water column that has N° = 0). 

Killworth (1980) showed that the character of the instability (if any) 1s 
determined by the behavior of the quantity 7670 /82(V,/N?). This quantity computed 
from the initial temperature field changes smoothly by 3 to 4 orders of magnitude from 
the sea surface to the bottom. Since the vertical scale for V 1s smaller than the vertical 
scale for T, this term approaches zero at great depth. This type of behavior for this 
term means that mixed (barotropic / baroclinic) instability 1s possible in the upper 
ocean lavers, but only barotropic instability 1s possible at great depth. 

We tried three different experiments to determine the optimal placing of the ten 
model sigma levels. Firstly, the original model sigma levels which included wind 
forcing, had finer resolution in the upper ocean than lower ocean. Secondly, ten equal 
spaced sigma levels which gives relatively poorer vertical resolution in the upper ocean 
were tried. Thirdly, we choose a distribution of levels which has good vertical 
resolution of 6/02(V,/N?) throughout the water column because the study focuses on 
the influences of topography. The final decision for the sigma levels is shown in Table 
2 


TABLE 2 
TABLE OF SIGMA LEVELS 


sigma Ol 035  O7S9s> Seer 
depth(m) 22.5 lOM38247 Sees 757.5 
sigma 315 441 597 ~~ .783 l. 
depth(m) 1192.5 1701.0 2335.5 3105.0 4011.7 





A plot of the cross section of the initial baroclinic velocity field used in the four 
experiments (computed from the eqn 2.1 using the thermal wind equation) is shown in 
Fig.2.7. The strongest velocity at the surface is 40 cm/sec (southward) at 40 km off- 
shore and zero velocity at depth 720m. This particular southward jet is representative 


of the observed CCS coastal jet in the spring and early summer, having been spun up 
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by upwelling-favorable winds (Kosro and Huver, 1986). Its transport, which is 5.1 Sv, 
is also representative of observations (Freitag and Halpern, 1981). The temperature 
field is shown in Fig. 2.8. The surface temperature is 7.5°C at the coast and 17.5°C 
off-shore which are also representative of summer observations. The potential vorticity 
is Shown in Fig. 2.9. The potential vorticity was computed from the model 
temperature and currents using the expression from Holton, 1981, page 92. In Fig. 2.9 
the potential vorticity is positive except near shore above 150m. The horizontal 
gradient of potential vorticity is shown in Fig. 2.10. Since the horizontal gradient 
changes sign in the domain, instability of the jet is possible. However,since the gradient 
of potential vorticity in the upper ocean is everywhere negative, the horizontal shear 


will probably serve as a stabilizing (B-type) effect. 
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Figure 2.8 Temperature field ( 
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Il. RESULTS AND ANALYSIS OF NUMERICAL EXPERIMENTS 


A. EXPERIMENT 1 (OFF-SHORE SLOPE ONLY) 

The first experiment considers the response of flow over an off-shore topographic 
slope which has the continental rise depth at the coast of 3000 m and the deepest 
ocean depth of 4500 m (Fig. 2.3). In this case the maximum surface velocity is 40 
cm/sec southward. The experiment of the flow over the topography neglects the depth 
average velocity component. Onlv the baroclinic flow over the topography is studied. 
The effect of this approximation will be examined in future studies. 

The experiment was integrated to 40 days to determine the response of the flow 
to the topography. The resulting surface U-field (eastward, or on-shore velocity 
component), developed in a time series (every two davs) 1s shown in Fig. 3.1 and Fig. 
3.2. The plot shows only the eastern half of the model domain, within 256 km of the 
eastern boundary. There are only small waves during the first ten days and they are 
entirely due to the initial random temperature disturbance. The disturbances propagate 
southward at about 14 cm/sec. The dominant wavelength increases with time from 
about 50 km at 18 days to 130 km at 30 days. There are only small amplitude waves 
north of 380 km. The eddies become stronger with time as they propagate southward 
(spatial and temporal growth). The eddy velocity decreases with increasing depth, and 
is very weak below 1000 m (not shown). During the entire experiment, the growing 
disturbances were tilted (in the (x,y) plane) in the direction of the mean jet, a tilt 
associated with barotropic damping of the eddies. 

Fig. 3.3 shows four different fields at 40 days. The in-shore temperature is 9°C 
and off-shore temperature is 18°C. The temperature field is affected by the geostrophic 
velocity which is determined by the pressure field. There is a strong baroclinic 
disturbance near y = 192 km, whose isotherm pattern has been distorted into a kind of 
“cold front”, propagating southward. The growth of this disturbance must be due to 
baroclinic instabilitv because, as noted above, the trough axis tilts in such a way that 
eddy V-momentum converges into the mean southward flow. The pressure field at 40 
days at depth 125 m (not shown) and surface are in phase, indicating that most of the 


baroclinic development is about ended. 
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A flat bottom case was integrated to 40 days to determine (by comparison with 
Experiment 1) the response of the flow to the topography alone. Fig. 3.4 shows the 
same four different fields for the flat bottom case at 40 days. By comparing the flat 
bottom case and Experiment | we can see that the flat bottom case has stronger eddy 
velocities than Experiment 1. In fact, the eddy velocities in Experiment 1 at day 32 
were about the same as those in the flat bottom case at day 22. Thus the continental 
rise topography has a stabilizing influence on the mean flow, so the flat bottom case 
goes unstable approximately ten days earlier than Experiment 1. As shown below 
(subsection E). the stabilizing effect is due to the topographic P-effect. This effect can 
be seen by the fact that the eddies in the flat bottom case do not tilt in the direction of 


the mean flow as much as the eddies in the topographic case. 
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Figure 3.1 U-field in Experiment | from 2-20 days in time series (every 2 days), 
with contour interval 2 cni/sec.. 


Dashed lines show negative values and indicate off-shore flow, while solid lines show 


positive values and indicate on-shore flow. 
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Figure 3.2 The same as Fig. 3.1 except 22-40 days. 
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Figure 3.3 Four different fields from Experiment | at 40 days. 
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Figure 3.4 Four different fields from the flat bottom case at 40 days. 
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B. EXPERIMENT 2 (ESCARPMENT SHALLOWEST TO THE NORTH) 

By comparison with the actual bottom topography (Fig. 2.2), the escarpment 
with the shallowest depth to the north (Fig. 2.4) is the most realistic, though still 
idealized, bottom topography for the California Current region. The width of the 
escarpnient is 60 km and the top of the escarpment is at y = 360 km. This means that 
the steepest part of the escarpment slope is located in the center of the domain at y = 
320 km (or latitude = 39.5¢ as in Fig. 2.4). The depth change due to the escarpment 
is 1500 m. The experiment was integrated to 40 days. The U-field time series (every 2 
days) is shown in Fig. 3.5 and Fig. 3.6. By comparing the surface U-field with 
Experiment 1, the patterns look similar during the first 10 days, but from 30 days 
onward a steady flow affected by topography, with symmetry at the center of the 
escarpment, is developed. The zero velocity is along the center of the escarpment (~y 
= 320 km) with off-shore flow (negative U) on the north side of the escarpment and 
on-shore flow (positive U) on the south side of the escarpment. The in-shore velocity is 
larger than the off-shore velocity at surface, but at the depth 790 m, the in-shore 
velocity and the off-shore velocity are the same. This development of a quasi-steady 
pattern, with no eddies, after 40 days is very different from the previous case 
(Experiment 1) having only off-shore variation in topography. 

Fig. 3.7 shows four different fields at 40 days. The temperature field at the coast 
is 9 oC, and offshore it is 18 °C. The temperature increases upstream (north) of the 
escarpment and downstream (south) of the escarpment with the lowest temperature 
directly over the escarpment. The surface pressure field shows the geostrophic flow 
pattern best. The pressure pattern shows a trough over the escarpment with a ridge 
downstream of it. The isobars are pinched together (as are the isotherms) on the 
downstream side of the trough indicating a strong jet-like flow toward shore. 

From the pressure field, it is possible to determine the effects of the escarpment 
topography with the flow moving from the shallow water column (on the north) to the 
deep water column (on the south). It is seen that the flow is stabilized by topography 
so Experiment 2 after 30 days, looks like steady flow. The topography induces a quasi- 
steady trough-ridge system over the escarpment. In addition, the time-variable eddies 
which were present in the pressure pattern of Experiment 1 (Fig. 3.3), no longer exist in 
this experiment with the escarpment. As shown below, the eddies are apparently 
stabilized by the strong topographic -effect which is strongest (largest topographic 


slopes) in this experiment. 
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Figure 3.6 The same as Figesresexecemenpenmenie:, 





P at day 40 deptho 





DUOCTOTOOEEnNG 
DLETEULIM EGG 
PEUEQUOUISIONNIET 
POU TDOEOMOOMOIEGLE 
CUOORUUTIEGET 
CPUC OUUOOGOT 
DOLD OCTET FEN EE 
TEIN 
4 
I 
f 
' 
' 
j 






COCUIMIMEE 
VEO TTEOHE GOON 
PE 
PE OVEEINTINAGSE 
STOUT OMEN 
Veueeanimins 


DUCT et 
Pee parerienimmesr 
setae fee rereeiwil 
of fee eee AR 
a pf of ee rereerenait 
tog ff fe eeeereaan 
Oeg ff fs taetrteonaien 
Seg J Ze tren 
' fsesveeryaperen 

ft rervsvaniy 

Or vovariynny 














a 
a 
> 


al 

ho 

oO 
Distance (km) 










43 






PeUsatasenane 
"Ute ripen 
VU detersnais 
pS LEeSerried 

Supa sanenuesas 


4 
SN S9ariassentii 
4 UNAS NNN ATES 
VEN Sa UNSEAT 
‘ ° Say Sh sttiisst 
ase setts 
uy Soy isons 
alti dann 
PML TT 
i ‘ Tiny 
RY Say EM anisttea 
\ tgtn eT 
Lats. "| mittay 
vy 4 














We. po Ht 
8h ts a@ MGT 


ch. . 
| i! vi roo O0NLy 









i! 127 Fey ame 
att ttny ‘even 
Witte (iter 
uty nyunranen’ 
ily thisanens 
' rs Uy, terete! 
toett Neeemaene! 








512-48 38k" 320 256 “182 "128" 0 512 448 384 320 256 192 
Distance (km) Distance (km) 
contour interval= 2.0 contour interval = 


T at day 40 deptho 























128 64 6 


1.0 


U at day 40 depthO V at day 40 depthO 


640 
' ' 
' 1 
' 1 
576 ' \ 
\ I 
I 
§12 : 
i 
I 
448 
os 
384 € 
<x 
— 
@ 
320 9 
‘= 
oO 
tet 
2 
256 5 





a 
512 448 384 320 256 192 128 £64 o 512 448 384 320 256 192 
Distance (km) Distance (km) 
contour interval= 2.0 contour interval = 


Df UPAD SPEC CUE OOO eH 
D DOUPEEE CE OUOOTIE OSE OTSt@N 


DEPEOUE EC EEEEIDECUOT TERI wha 
DUCUE EDAD ODDEST Ebi ete war 
DIDODEEE OTE Ie wa 
DOV TUEO UCP EEE SOME tteeD ia 
DULCE TEED ETEEDO MD OnNetenes of 
PUP EE OU EP UOTE DOGO SET tenes Fd 
POPUP EDC OEEDODT Ee Ya) 
PAPEESULECAOEEOSEELaT em OARS 
POOCCUU UO DODEOOOIEEDAE mnie 4 
ml iaA 


} 
un Ce) 





UUCGNGANOT DING WAIT 
SN GSUGQQTTEISSINWIHID 
NV VSDUUyeaaapersiavadt 
NT VAG Arg gagsusingiig, 

VUSTSUUNDISSNGIORmIN 


LAPRP CVU DD 


ERURIELTSTUC yy 
Sof POIs Tatas 
PP UUTP TN CASTS y 2 
SOOT UL LNLUMEERINY RG 
4s PeL IIA 
Nesey tO 


= 


ar 


Sete 


- wT 
pean feet 


Ve ace 7eeeeei 
pless 


‘1 hd 


29 ybC gM! 
Hae 
tf Vea) 



















—— 
38a c 
x6 
— 
oO 
320 oO 
Cc 
o 
“ 
256 = 
192 
128 
64 
a 
640 


576 


512 


448 


ia 
oo 
> 


Distance (km) 


320 


256 


192 


128 64 6 


2.0 


Figure 3.7 Four different fields from Experiment 2 at 40 davs. 



























C. EXPERIMENT 3 (DEEPEST TO THE NORTH) 

The escarpment of Experiment 3 (deepest depths to the north) is the reverse of 
Experiment 2 (shallowest depths to the north). The top of the escarpment is at y = 
280km (maximum slope at y = 320 km) and the width of the escarpment and the 
depth change due to the escarpment are the same as in Experiment 2. This experiment 
was integrated to day 40. The surface U-field time series (every 2 days) is shown in Fig. 
3.8 and Fig. 3.9. The eddies become stronger with time as they propagate southward. 
There is a large off-shore flow at about y = 280 km. The in-shore flow and off-shore 
flow is not symmetric across the center of the escarpment.There are only small waves 
at the north of the escarpment and these waves grow rapidly at the south of the 
escarpment. 

The surface temperature range (not shown) is the same as in Experiment 2. The 
temperature field is affected by pressure which indicates the geostrophic velocity. The 
southward surface flow moves from the deep water column to the shallow water 
column. Particles displaced into shallow water tend to turn to the right in order to 
acquire anticyclonic vorticity, thus conserving potential vorticity. Although there is a 
strong tendency for the surface flow to be directed off-shore (U less than zero) over the 
escarpment, the eddy field is still quite strong at 40 days and the flow is not yet 
approaching steady state. The results of this experiment are very different from the 
case with the shallowest depths to the north. Thus, the orientation of the escarpment, 
whether the shallow water lies to the right or left (looking off-shore), makes a big 


difference in the flow pattern. 
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Figure 3.9 The same as Fig. 3.2.except Experimenta: 
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D. ENPERIMENT 4 (RIDGE) 

Experiment 4 uses a combination of the topographies of Experiment 2 and 
Experiment 3 to form a topographic ridge (Fig. 2.6). The top of the ridge is in the 
center of the domain at y = 320 km. The half width of the ridge (Gaussian length 
scale) is 60 km, so each side of the ridge is exactly the same as the similarlv facing 
escarpment in the previous two experiments. The U-field pattern (Fig. 3.10 and Fig. 
3.11) looks like Experiment 2 but the eddies are stronger than Experiment 2. The 
eddies move to the south and get stronger with time. The in-shore flow downstream 
from the topographic ridge is larger than the off-shore flow upstream of the 
topographic ridge. The off-shore flow is over the center of the ridge. This flow is also 
consistant with the conservation of potential vorticity. Potential vorticity would 
increase when depth decreases as particles approach the ridge. However particles turn 
to the right (off-shore) and acquire anticyclonic relative vorticity thus maintaining 
constant potential vorticitv. The reverse happens (particles turn to the left) as the flow 
goes over the ridge. The result is a trough (positive relative vorticity) over and 
somewhat downstream of the topographic ridge. From the surface U-field, we can see 
that this experiment looks like quasi-steady flow after approximately day 30. This 
result is very similar to the escarpment (only) case with the shallowest water to the 
north (Experiment 2) but very different from the case with the shallowest water to the 
south (Experiment 3). It therefore seems that, if the variable bottom topography has a 
southward facing slope (westward topographic p-effect) as in Experiment 2 and 4, a 
quasi-steady flow pattern develops. On the other hand, if the topographic B-effect has 
no westward component (Experiment | and 3 and the flat bottom case) then no steady 
state is reached within 40 days. The four topographic experiments are compared and 


interpreted in greater detail in the next section. 
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Figure 3.10 The same as Fig. 3.1 except Experiment 4. 
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E. ANALYSIS AND PHYSICAL INTERPRETATION 

The above results can be seen to be characterized by two main features - (1) the 
establishment (or not) of a quasi-steady response to the topography after 40 days and 
(2) the growth (or not) of the initial small amplitude random disturbances. In this 
section we try to interpret these main features. 

The four different experiments with variable bottom topography were integrated 
to day 40 and the four resulting surface pressure fields developed in time series (every 
ten days) are compared in Figs. 3.12-3.15. The four different fields look much the same 
at ten days (Fig. 3.12). By comparing the pressure field pattern at day 40, Experiment 
3 and Experiment | look similar except the amplitude of the main eddy in Experiment 
3 {the escarpment deepest to the north) is larger than in Experiment | (off-shore slope 
only). Also, the surface pressure field in Experiment 2 and 4 look the same. The 
maximum off-shore velocity in Experiment 2 is at the top of the escarpment (y = 360 
km) and the maximum off-shore velocity in Experiment 4 is at the top of the ridge (y 
= 320 km). Experiments 2 and 4 have reached a quasi-steady state but Experiments | 
and 3 have not. 

Some of these differences and similarities may be explained in terms of the 
direction of eddy energy propagation. Due to the topographic f-effect, the eddy energy 
of the quasi-geostophic, non-dispersive long topographic Rossby waves propagates 
with the shallow water to the right. The eddy energy in Experiment | therefore 
propagates to the north and that in Experiment 3 propagates along the escarpment in 
the in-shore direction before turning northward near the coastline. The eddy energy in 
Experiment 2 propagates along the southward facing escarpment in the off-shore 
direction, and that in Experiment 4 propagates in the in-shore direction on the north 
side of the ridge and in the off-shore direction on the south side of the ridge. Thus, 
there is a mechanism which tends to converge the eddy energy toward the eastern 
boundary stronger in Experiment 3 than in any of the other experiments. The eddy 
energy can therefore propagate away, leaving the steady (topographically forced) flow 
behind, most easily in Experiment 2 and 4 because in both of these cases, and only in 
these two cases, there is a westward off-shore component to the group velocity vector. 
This is why the orientation of an escarpment, whether the shallow water is to the right 
or whether it is on the left when looking off-shore, is important. 

We now compare the growth of the initial random disturbances in the different 


experiments. The flow in Experiment | acquires an eddy velocity of 30 cm/sec at about 
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day 32 and that in Experiment 3 does so at about day 30. The flat bottom case has 
eddy velocities this big (at about day 22) earlier than Experiment | and 3. These 
differences in the growth of the random disturbances are probably due to the different 
(stabilizing) effect of the topographic-B in the experiments. We estimate the effects 
using the simple two-layer model. From Holton, (1981, pages 216-223, two-layer 
model), we use two parameters B and A to compute the phase speed and growth rate of 
the disturbances. The parameter B/h is the mean potential vorticity gradient on an f- 
plane in the case of a small Rossby number (i.e., topographic-B), and d is the inverse 
of the internal Rossby radius of deformation. These two parameters are defind as 


follows: 
M? = Poi(N*Oz*) ~ 2.0 E-13 fem? 
B = fh X Gh/6x ~ 2.5 E-12 1/(cmsec) (for experiment 1) 
B ~ 5.0 E-12 1/(cm.sec) (for experiment 3) 
V_, = 20 cm/sec the vertical average meridional wind (upper 720m see Fig.2.7) 


V_ = 20 cm/sec the basic state thermal wind (from Fig.2.7) 


Pele) kin meridional wavelength of eddies 
k = 2U, ie = 5.23 E-7 meridional wave number 


aie B(k> + 27)/k7(k? + 24°) & Jd phase speed 

From Experiment | (continental rise) we obtained a phase speed (C. in Holton) 
of 11 cm/sec southward and an e-folding time (1/kKC. in Holton) of 2.8 days, and for 
Experiment 3 we obtained a phase of speed of 2 cm/sec and e-folding time of 3.4 days. 
From the model results in this study, we estimated the phase speed and growth rate of 
the eddies from the resulting U-field (Fig.3.1). We obtained a phase speed of 14 cm/sec 
and growth rate of about 10 days for the Experiment 1 and about the same for 
Experiment 3. Thus, the random disturbances are baroclinically unstable and grow at 
about the same rate in both Experiment | and 3. By 40 days, the flow in Experiment 3 
is in the very early stages of the approach toward a steady state flow over the 
escarpment (on-shore / off-shore flow). The disturbances appear to be somewhat 
stronger in Experiment 3 than in Experiment 1, not because of a stronger instability 


(shorter e-folding time), but rather because of the convergence of Rossby wave eddy 
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energy at the coast. In Experiments 2 and 4, the topographic B-effect appears to be 
strong enough, and the direction of eddy energy propagation away from the coast, so 
these two cases have already reached steady flow by 40 days. The resulting new “mean 
flow”, ie. the quasi-steady trough / ridge system over and downstream of the 
escarpment (or ridge) in these experiments, appears to be stable. Future experiments 
should be aimed at examining the stability of this topographically forced flow to a 


variety of disturbing effects such as variable winds or boundary forcing. 
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Figure 3.13 Four different pressure fields at day 20. 


640 
576 
Bre 
448 
— 
384 £ 
x 
ee” 
o 
320 UY 
= 
Oo 
D 
256 5 


192 


128 


64 


640 


576 


312 


448 


Gs 
@ 
> 


tA 

BRO 

o 
Distance (km) 


KD 
an 
a 


192 


128 


64 


44 






Distance (km) 


contour interval = 


SPI RII 
00705 LA eam | 
SRA ih 
SO gy OFM f 
(ties opt artpatnmepeyl 
FONTS Sime, arn 
Vl 09000000, 7° 4g 
SVaderee 
SSsahont 


512 448 
Distance (km) 


contour interval = 


P at day 20 deptho 


TEC APL SAT M28 @@ mei) 
LULU Ct aeees wel) 
JIL BHetem ee ee) 
(PCI DINee om ee er 
SLAUU PULL ee ee ma ia 
#122 LNA wewrneney 
FU Clepesed ae an) 
eT Nee eae 
feet Hillman aaeneti! 


einsa™? e'g! efi ». 


og tp thel ov 

Weseted wt gine 
V UL ETT TT Tap ty tee oll 
SUIC CF LIN aa atrg ae 
SSSVIOL CIN OReOmen) 
SULUIDD @2e7 RRR pore 
CU INULL Wee 2 srt bE 1 1 
UOC PLLL IN tt4/72 eat 
VPP a see as nee; 
Flee enttreaea|, 
OEFCLEOEEE LOMA AOI, 15 | 


4 Orta a ts ENO Etat; 
PrP esyilel lie Rea 


PE soca nilMateng, | 
OFF bereeeesirotinng 


SU UG ceepeneeee ee aiatery 
mil 
PVUN a peer leeeeer Aan gage 
SSS senate IIe teatag t 


S88 sop yyyeel eh cel LAD naa) 
. Hin 
{Pawan AY 
\wauneas! 


Als 


[Cine 
SConyeeses: : 
C1yga atl ‘tae 





' 

' 

ligt ta 
{ eT eer) 1 
ty Ee aa caneenn en 
1 


1114 atte 


128 64 0 





iO) 


at day 20 deptho 


HUIS INI peemeee 
CSTS0TIIG 0a me) 


LUSULIIIIMI0@ commen 
PS ULTIO0S)0 @ Nao ait 
Distliinhe heap 
LUSSS1S1 Astor: 
VS 0001) 1a 
C8 08800 preg 32) 
109990 900 Pirates 9) 
LUV1 1800S 
40000 901 (o cp 
18009000 


a 4 





' 
‘ 
Vas 
vil 
ile 


éline 


I 
' 
' 
tt. 
iv 
i 
@lelfwee 

UD 


(Ohi, 
aff stews 


Pe Ciiwe@eaum 
a} I Ar, 
SCT or 

i 


“ 


ws, 





no 


640 


576 


312 


448 


da di 
N a 
o - 


PO 
in 
mH 
Distance (km) 


640 


376 


312 


448 


tA Gi 
ne oo 
o » 


AR 
in 
aH 
Distance (km) 


192 


P at day 30 deptho 


640 
CUO FOO a 
CUO UCC CeNeCre ster net 
VU CU OLN eeeeee cmon 
VU CU SCNCC eee mCeee 
ECCOFOT ANNO 
OCCU CC ernO tT 
CO UCC OOM onennenenetoreNE t 
COSTE CONSE NEmeottALe 
VOU NN eanceeeettsnamengt 
CCU N Dennen concatenate tet 
COU ONe TON se 

CU UON POC Ct Rae S22 
POUR 
PSU UON CNR Nao 
o'lelaliatsta! 
$2 Py Meta U yt 
Pdoceds WE yp ys 
EDERAL att M eae b 

PUPS TL eee 

Ved ts ee ceeeereranintay 

ULE De ateteesennmnnnsy 4 

FUVNL OP Vormrnte sey | 

ANETTA LUMO e Ze 
SATIRE E To ial 
Savanna sreense Onan 2 
NVUU Supe ges sections | 4 
UVM ti UMAAMUV A. NS 
SUVS aeuvupneemmmeaays S 4 
UU taney vwcmenee yy 
DTV EATT EY Ub aah 
SVevtec pee won 
CPTI 1 doptiedtttedlL 
448 £1ayttthaeeAnsss 
O78 Se 2; ggg; 
PATE O OS, gg AMEN Y 

6erebtran gy Z 

£12 ed arene ers 

ee OU ES 

SAGE EES prema RUAN 
VANE ESA yen AB SS 
PAN ALEVE N a 


. 
UA AERA ul 


AN i 
LASS aNag geese’ 
f : Wahstya ae can tee, 192 


576 


448 


G 
a 
- 


Ga 
N 
oO 


Le] 
in 
® 


128 


a 

‘ # htuecaeet 64 
2 aanneee “iy @ 

? 

; Pr nelrns, at 
tenes vi 
sol sw gee 
tal ial 4 


320256 mI92) 128" "64°." O 
Distance (km) 





512 448 384 


contour interval= 1.0 


P at doy 30 deptho 
640 


CPUC OURO E 

COUT ON eedenrsearwtereed tb 

LL eet ddan tectrreae mut ¢ 

PL pee ties musmeuats 

LUD Anon seat eec@usmeen tt 

VOC intiunmuie 8376 

DUCE PONHINOOMU ARMIN CE 

VAT NNDUUUONIOOEOAHURIOD 

POL t ene 

CE ROOND OI emumerbaS 

1 LCL ON StrteAneett 
HLiMiiMeu gO 12 
VAGOMsiCMAeseUts 

DEAR AY 


’ 448 
Rd) 
07 0) voemamiames 
V0: (Ae 
AVecd ii tyeeetmmy 
NAV SUD ND eeeteseey 
SUAGUU UND IAM RR ¢ 
SVOON Ce Up OO mma; 
VConiiectaasen, 
Ch onnie eee, 
PA UUme 4) 
lia taurine 5) 
22a, ge pene ana) | 
Hh ATH | 
Opt Zl HRN | 
PS Oe FTN | 
tetsin f 640 ea § 7 
tine, « 
yal)? 
snt,? 
' 
vse 


vaaatslee 


' 
' 
' 
4 
4 
é 
' 
‘ 


Ga Ga 
nN oo 
° a 
Distance (km) 


NO 
an 
m 


é 
{CAM pepeeey 
ws 
‘ wt Hams ~~, 
4 wi Wee 


leyh y a 
Ny ap mreren . wey 192 


, 
Netti wens 
RSL ae 

SQV aan’ yh. 12 8 


‘ \ 
is Mamvwunmiusl® > 
wee anes 1 
ys runs 
x Sy x 
y 8 2 
\ 
wh 





512 448 384 320 256 192 128 64 o 


Distance (km) 


contour interval = 1.0 


Distance (km) 


P at day 30 deptho 


S12 448 384 320 256 192 


Distance (km) 


contour interval = 





640 


RTT 
Chenin 
(Mtns 
(1 ¢te 00001 oT al 
CeCttCte00) a oo 0: 
COCtttt00c 0 oo 8 
C1 UCC C0010 60 90 eo 0 
Me 
JU010000100m oo | 
104 C01 1000000 eee 
rwiiinema:_ =r 5 12 
#0 10000(200010 6 mmm 

211000240 e mE 





576 


e fy oO fe tg et alow! 


sb ty teal te 
Pred a i an 

O00 M Ot eet) 

ARP 
PACs He 
MINH Amen) 
POUR WAAROE) 
lelisaseAdAleen 
(410 Cert eiaemmen 
(P0000 e002 LHe! 
Uren aren 
PON eet AL ew 
FAVE tre Arel 
POPU AURA 
Pelee seeuRiinTimena 
SVFNCC Cleese MAA 
SVAN AMO e@ meee) 


448 


ee) 
a 
- 


4¢ 
ee 
ue 


A 
No 
oO 


va 
4 


PQ 
cn 
ni 
Distance (km) 


192 


128 64 ® 


1.0 


P at day 30 deptho 





S12 448 


384 320 256 192 


Distance (km) 


contour interval = 


640 
19000140) cea 
1IS@0(0( 
CU C00 (\0(}.o 
1 
(17 0100Ve Oo 
FINN 
cero 

PV GUUihe 

Cf test 

cai 

VUDC1\ (foe 
NAc a & 12 

Pecuiten 

cunsene 


576 


448 
tf euing 
yj 
ff licma 
ee utunu 

deta 

fei 

4 eee 


G 7) 
N a 
Qa eS 


RQ 
tn 
m 
Distance (km) 


192 


1.0 


Figure 3.14 Four different pressure fields at day 30. 
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Figure 3.15 Four different pressure fields at day 40. 
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F. SUMMARY AND CONCLUSION 

The objective of this study was to investigate the topographic effects on coastal 
jet instability. By comparing Experiment 1 and the flat bottom case, we can see the 
continental rise topography has a stabilizing influence on the mean flow, so the flat 
bottom case goes unstable approximately ten days earlier than Experiment 1. If the 
variable bottom topography has a southward facing slope (westward topographic - 
effect) as in Experiment 2 and 4, a quasi-steady flow pattern develops. If the 
topographic p-effect has no westward component as in Experiment 1 and 3, then no 
steady state 1s reached within 40 days. From a careful and complete analysis of the 
analytical and numerical experiments describe above, we conclude that Mendocino 
Escarpment plays an important role in baroclinic and barotropic instability processes in 
the: CCS. 

Since the depth averaged flow was not included in the model used in these 
experiments, the above studies should be repeated using a model which includes the 
depth averaged flow. This model could allow the above conclusions to be reevaluated 
and it would reveal the role of the depth averaged flow, and its interaction with the 


baroclinic flow and topography in the CCS. 
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